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Abstract 

In this paper we put forward a mechanism in which imploding shock waves emit electromagnetic 
radiation in the spectral region Ao = 2itRo-, where Ro is the radius of the shock by the time it is first 
formed. The mechanism relies on three different pieces of Physics: Maxwell's equations, the existence 
of corrugation instabilities of imploding shock waves and, last but not least, the Inertial Polarization 
Principle. The principle is extensively discussed: how it emerges from very elementary physics and 
finds experimental support in shock waves propagating in water. The spectrum of the emitted light 
is obtained and depends upon two free parameters, the amplitude of the instabilities and the cut-off 
i?max, the shocks' spatial extension. The spectral intensity is determined by the former , but its 
shape turns out to have only a mild dependence on the latter, in the region of physical interest. 
The matching with the observed spectrum requires a fine tuning of the perturbation amplitude 
£ ~ 10~^* , indicating a quantum mechanical origin. Indeed, we support this conjecture with an order 
of magnitude estimative. The Inertial Polarization Principle clues the resolution of the noble gas 
puzzle in SL. 
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The Inertial Polarization Principle 



In this paper we put forward a mechanism responsible for transducing the kinetic energy stored in an 
impfoding spherical shock wave into electromagnetic radiation, which is based solely upon Maxwell's 
equations, the existence of very small instabilities away from the spherically symmetric flow and the 
inertial polarization paradigm. Based on these premisses we obtained the spectral intensity of the 
outgoing radiation. The mechanism turns out to be so efficient that the observed energy emission rate of 
P{\) ~ 10^ ^'^Watt/nm calls for perturbation amplitudes no larger than e — 10~^^! Maxwell's equations 
are a pillar of theoretical physics while inertial polarization is a consequence of very elementary physics 
: an atom that undergoes an acceleration, say a , develops in its interior polarized electromagnetic 
fields. The issue is made clear for an observer sitting in the frame of the molecule, where he sees inertial 
forces acting both upon the nucleus F/v = Mj^a and on the electronic cloud Fg. = M^a. The gradient 
between these forces tends to sag the cloud away from the nucleus, and the atom develops internal 
polarization fields, say Eq, to compensate this gradient e Eq ^ (Mjv — Me)a . The role of inertial 
polarization remained hitherto unnoticed only because detectable polarization fields call for tremendous 
accelerations, say, Eq ~ IV/m would require a ~ {e/Mp)Eo ~ 10'"^ {eV/{MpC^)){c^ /cm) ~ 10^"cm/sec^ 
which are absent in every day life experiments. Nevertheless, there are two instances where such large 
accelerations manifest: i.) in the realm of very strong gravitational fields where inertial polarization was 
shown to be the working mechanism that rescues the second law of thermodynamics from bankruptcy 
(otherwise super-luminal motion of black-holes inside dielectric media would entail a violation of the 
generalized second law ii) in the realm of shock waves, because shocks are powerful accelerators 

of fluid molecules: a fluid molecule that crosses the shock undergoes a macroscopic velocity change (of 
the order of the fluid velocity itself ) within a microscopic distance - the shock width (of the order of the 
mean free path for the atomic collisions |^). 

The inertial polarization principle is the single non- very-well-established piece of physics in our recipe 
and we proceed by making our case for it. Consider a planar strong shock wave propagating within a 
perfect gas. Let V2 and vi represent the fluid velocity in the back and in front the shock, respectively 
( likewise, the index 2 (1) refer to physical quantities behind (in front) the shock ). As the fluid molecules 
cross the shock they experience a mean acceleration a — {v2 — wi)(A<), where At is the mean time it 
takes the gas to cross the shock-width S. Clearly At — 6/v, where v ^ {vi + V2)/2, is the mean velocity. 
Putting these pieces together 
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For a strong shock propagating in a perfect gas ||] : 
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where Vi is the gas' specific volume. The compression rate satisfies V1/V2 = (7 + l)/(7 — 1)[| 
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The shock width 5 is known to be of order of the mean free path for collisions of atoms in the fluid, 
S ~ (na)~^, where n stands for the number density of atoms and a for the collision's cross section. 
Bearing in mind that nV = A/ ^ where A is Avogadro's number and is the molecular weight of the gas, 
we obtain the colossal figure for the mean acceleration atoms experience as they cross the shock: 
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The mean electric polarization developed across the shock Eq « {Mp/e)a is also sizeable 
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Unfortunately, the shock is so thin that the voltage developed across its ends is very small 

1/ . ^ ^ 1.2 X (J^) ( (6) 

e 7 — 1 V atmJ \ g Q J 

Shock Polarization was first observed in the early sixties Q for shock waves propagating inside water. 
Since then, both quality and range of the measurements improved considerably j^. Harris |^ credits 
the effect to the fact that large pressure gradients inside the shock results in a torque field acting upon 
the water molecule causing the molecule's dipole to align. We reproduce his results via the table: 
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The underlying Physics for a shock propagating in water is the very same as for a gas and we infer 
the averaged electric potential across the shock from eq.@ bearing in mind that: i.)~the compression 
rate for water is of order one, therefore we go one step back in this equation by replacing 7/(7 — 1) — *■ 
7/(7 + 1) ~ 1/2); ii.)"'the equation was obtained for a gas and for liquids it should be regarded as the 
linear expansion of the function V{p). Then it follows that V{mV)/p{kbar) ~ 0.6 , in agreement with 
the lower pressure region of the experimental data. Detection of shock polarization for non-polar fluids 
would vindicate the Inertial Polarization Principle. 

The acceleration field inside planar shocks is space and time independent (and so the corresponding 
polarized electromagnetic fields) . Nevertheless, planar shocks are known to develop corrugation instabili- 
ties ^ , small deformations of the planar geometry that detach from the shock and propagate throughout 
the fluid. They correspond to the spontaneous emission of sound from the shock. These instabilities will 
cause a space time dependent acceleration field inside the shock, and by the Inertial Polarization Principle 
a wiggling 6 x lO'^V/m electric field vector that is radiated away: sound and light are emitted simulta- 
neously , provided the Inertial Polarization relaxation time is small enough. This brings to one's mind 
the famous and intriguing sonoluminescence effect 0| in which under heavy bombarding of ultra-sound 
waves, a little ( 5/zm) bubble of air cavitating within a flask of water undergoes a spectacular collapse, 
attains the supersonic regime and glows (mainly) violet light. The effect has been around for sixty years 
or so ( @,[?]) and proper understanding of the problem remains elusive. The most popular mechanism 
is the Bremsstrahlung from free electrons in the gas where the ionization is caused by two successive 
heating processes: flrst the adiabatic collapse of the bubble which is then followed by the motion of a 
shock wall inside the bubble (the shock's Mach number controls the temperature rate T2/T1 ~ M^) pof . 

The formation of a shock wall, a collapsing spherical front of radius R{t) = A{—t)°' (a < 1), happens 
by the time supersonic regime is attained inside the bubble Q . The acceleration of the shock front surface 
a{t) ~ A{—t)"~^, becomes very large at focusing {t 0) engendering very large space and time depen- 
dent Inertial-Polarization fields. Nevertheless, the spherical symmetric geometry of the problem prevents 
these fields to be radiated away: pursuing the present avenue seems to require some supplementary mech- 
anism to account for the radiation flash (a sparking mechanism was proposed |jll|, |l^). Fortunately, no 
supplementary mechanism is needed: numerical calculations (p^) have shown the existence of unstable 
perturbations of the collapsing shock which provide the multipole time-dependent inertial-polarization 
fields that are radiated away. The purpose of this paper is to calculate the spectral distribution of the 
emitted light . 

The paper is organized as follows. The following section reviews the dynamics of imploding shocks, 
and the existence of unstable multipole perturbation modes is rigorously proved. As a bonus, we obtain 
the energy and the power carried away by the sound waves that detach from the shock (corrugation 
instabilities). A novel semi-analytical procedure for solving the differential equations for the perturbations 
is developed, which nevertheless, is displayed in the appendix in order prevent the disruption of the 
main argument line with technicalities. In section II , we obtain the polarization fields engendered by 
the corrugation instabilities and show that they act as a source term in Maxwell's equations. Then 
we calculate the spectrum of the outgoing radiation. The spectrum depends on the dynamics of the 
corrugation instabilities, but fortunately it is possible to obtain the main structure of the spectrum 
without having to delve too deeply into the dynamics. The intensity of the outgoing radiation turned 
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out to be proportional to pe^ where e is the corrugation instability amplitude and p = Ep/{2ah), is 
Inertial-Polarization power-constant {Ep stands for the proton's rest energy and a for the fine structure 
constant). This constant is of the order 1.47 x 10^^ Watt (!): collapsing shock waves are the most 
efficient power-stations in nature, with the sole possible exception of astrophysical objects! Agreement 
with the experimental data calls for amplitudes of the order e ~ 10~^^ or 6r ~ 10~^^m! These tiny 
perturbations must have a quantum mechanical origin, and we support this conjecture by an order of 
magnitude estimative. Finally we suggest the resolution of the noble gas puzzle in SL. 



1 Dynamics of Imploding Shocks 

The non-viscous implosion of a spherical shock cannot be characterized by any dimensional parameter . 
Consequently the flow admits a self-similar symmetry. Let R{t) = Ai{—t)" represent the radius of the 
shock front, where Ai and a are two constants and Vshock = aR{t)/t, its implosion velocity. The self 
similar parameter here is ^ = r/R{t); the surface of the shock is given by ^ = 1. Self-similarity constrains 
the form of the speed of sound, radial flow velocity and density p3 : 
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V2 = {^)V{0 (8) 

P2 - PoG(0 (9) 

When expressed in terms of the self similar quantities Z, V and G, the boundary conditions for a strong 
shock n • Vshock >> c read. 



The equations that govern the flow are the entropy and mass conservation laws and Euler's equation . 
They provide a set of non-linear coupled equations for G(^), V{^) and Z(^), which when solved for Z{V) 
and S,{V), yield the pair of equations 0] 
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where k = 2(1 — a)/(a7). Inspection of these equations reveals the existence of a singular point at 
Z = {1 — VY ( dV/d£, 00?). Clearly, all physical quantities, and their derivatives must be finite across 
the singular point, meaning that the conditions {3V — k)Z — V{l — V){l/a — V) — Q and Z — {1 — V)'^ are 
simultaneous to each other at this point, such as to keep their ratio finite. Call T4(a), Zc{a) the solution 
of this pair of algebraic equations. The parameter a is obtained by numerically integrating Z(V) from 
V = V{1) to Vc for different values of a until the matching Z{Vc{a)) = ZjyCx) is obtained. The good 
values for a are 0.688376/0.71717 for a monatomic/diatomic gas. The limit i — » 0_ corresponds to the 
shock's focusing time, after which the shock reflects and reexpands. For latter reference, we mention the 
asymptotic behavior V ^"^1°" as ^ ^ 00 

We are seeking now perturbations away from this flow. Let 8 = Sp/p be the contrast function and Sv 
the velocity fluctuation. The latter can be decomposed into its normal and perpendicular components 
6vn ^ n ■ V , Sv± = Sv — Svnft. 

The linearized mass and entropy conservation equations read 

^+v^)s + 6v,.'-^ + V.Sv^0 (13) 
at or J or 
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while perturbing Euler's equation yields 



d d 



^dt dr 

Next, we introduce the self-similar ansatz 
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where Cp is the specific heat of the gas , Iq is shock formation time and e the amplitude of the perturbation 
at this moment. After some tedious algebra we translate the previous equations in terms of the self-similar 
quantities. The mass and entropy conservation yield (p^[?]) 

(l-F)C(A'-$') =/3A + 3$-/(? + l)r (17) 

(1 - V)ia' ^ I3a-n^ (18) 
where k — 2(1 — a)/ (07). The projection of Euler's equation (^5|) into the perpendicular direction yields 

1, 



a compact form 

(1 - V)^t' ^{2V + (3- -)t + Z{A + a) 
a 

but the normal projection gives a more cumbersome expression 

{1~V)^{{1-V)H' -Z{A' + a')) 
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(20) 
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Equation ( [l7| ) suggests the definition of a new dynamical variable 11 = A — $. We display these equations 
in matrix form 

±\Y{V))=M{V)\Y{V)) ; 0<V<V{1)^V^ (21) 
where \X{V)) = {(t){V),T{V),n{V),a{V)), ; \X{V)) = exp[/3 J^^ m{V)dV]\Y{V)) and, furthermore 



M{V) m{V) 
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■P[V) 
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03 (v^) 



Z[5 - 2/a -f 2/3 + (7 - \){'iV + dV)] + (1 - Vf[-l/a + 2V + 2dV] 
-Z + 

Z[{-f - l){iV + dV - k) ^ P) 
Z[j{3V + dV-K)+p] 



(22) 

(23) 
(24) 
(25) 



where dV{V) = . Clearly this set of differential equations possess a regular singular point when 
Z - (l-V)^ ^ 0, that is to say, at K. The limit V ^ ^ oo , m{V) -a/V, P{V) -> 1; 0i 
— 1/a; 1/12.3,4 0), reveals an additional singularity 

/ 1 \ 

-0(3 + alil\l) I I 
\ aK 0/ 

The matrix on the right-hand-side of this equation defines an eigenvalue problem whose solution 

\0i >= (0,0,1,0) 
\02 >= (0,0,0,1) 



Ai,2 = 



_ ( \e3>=il{l + l),3 + f3,0,aKlil + l)) 

- ^ { |04>-(o,i,«;(? + i),o) 

yields the asymptotic form 

\X{V)) ^V-"^[iai\ei> +a2\02>) + V{a3\e3> +ai\9^>)] ; V ^ 0, (28) 

where a„ are integration constants. Asymptotically regular fields require K(/3) < (except for the the 
particular mode ai = 02 = which calls for a less stringent condition 5ft(/3) < 1/a ). A further constraint 
on /3 arises from energetic considerations. The energy of a polytropic gas is 

E = / p[v'' + -\dV. (29) 



7(7 - 1) 

The lowest order contribution ( in the perturbation parameter e ) to the energy stored in the perturbed- 
shock is the second order expression 

5Ei [t) ^ 1 5p[v5v^ + Unr'^dr, (30) 

J 7(7-1) 

or after some algebra 

5Ei{t) ^ 4^a2eVoi?^^^(5^^Q, (31) 

where 

ci = G{am + m)W{i - v)^ + + h- imo + mm'dc m) 
Ji 7(7-1) 

and Rq stands for the radius of the shock by the time it is first formed to. Note that for ^ >> "1 , $(^) + 
n(0 V-""^ ~ C'^,G(0 const: the integral diver ges as ^5-1-2/3^ vindicating the introduction of the 
cut off ^c, which represents the boundary of the self-similarity solution. Clearly, this energy has to 
remain finite at any time and at focusing it requires that 1/a — 2.5 < Re(/3) < 0. In the appendix we 
develop a semi-analytical method for solving eq.(pTl) and obtaining the correspondent spectrum for „. 
In consonance with previous numerical calculations (p^) we confirm that (3 lies in this interval. By the 
way, the most unstable modes are shown to lie in the interval. 5 + 1/a < Re(/?) < —2.5 -I- 3/(2a) , even 
for very large values of /. For these modes, the energy emission rate 

i2a/3+5a-3 

Pi{t) = 4^(2a/3 + 5a- 2)a^e^p^Rl ^^^^^^^ Ci (33) 

diverges. This means that, in analogy with the corrugation instabilities in planar shocks, a burst of sound 
is emitted at the focusing. The total energy carried away during the shock-collapse is 

Esoun, = E ^mn) = (34) 

where we defined C = Re[5^;^-^ ^ C/ (/?)]. 
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2 Inert ial Polarization At Work 



As discussed already, electromagnetic bounded systems whose constituents have sizeable mass differences, 
say AM , and which are subjected to a strong acceleration field dv/dt engender polarization fields Eq, Bq 
that tend to restore the balance between electromagnetic and inertial forces. Clearly, these polarization 
fields satisfy 



dv 

AM— = Ze[EQ + -xBQ 



(35) 



where A and Z correspond to the atomic and proton numbers and e is the electronic charge. Clearly, 
AM « AMp, where Mp is the proton mass. Defining a polarized potential- vector ($0)^o) in the usual 
way, allows us to write the balance equation in the form 



--.x(Vx.) + V- 



Ze 
AMpC 



OAq 
dt 



-vx{V xAq) + V(c$o) 



(36) 



that suggests the identification Aq — v and <J>o — ^^J" i?" jl . Other possible identifications 
exist, but they are gauge equivalent. The corresponding polarization fields are 



_AMp 



dv ^v^ 



; .Bo 



— X V 
Ze 



(37) 



The time varying inertial-polarization fields engender the radiation fields i?, B and their superposition 
must satisfy the sourceless Maxwell's equations: 



V • (^; + E^) 

y-{B + Bo) 
1 d 



Vx{E + Eo) + -^(B + Bo) = 

W x{B + Bo)~-^(E + Eo)=0 
cot 



V ■E = AlTQejf 

V-B = 

c at 

IdE 47r,-> ^ , 

VXB--— = -(Je,/+ Je//) 



(38) 



with 



Qeff = 

Teff = 



AMp 
'AttZc 



dV -v 

dt 



AMp 
AnZe 



d^v 1 d ^: 



clearly satisfying the conservation equation dg^f f/dt+^- j ^ff = and 

J = — X V X V 

AireZ 

For non-relativistic flows ~ {L/T)'^ /c? ~ v"^ /c^ , and the field equations reduce to 



V xE + 



V -E 

V-B 
IdB 
c~dt 



VxB--— 
c at 



(39) 



= 



r/V X V X fT 
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where rj = AMpc/Ze. Next we expand ^ ^ J = y j where e„ is the familiar vector basis 

[|6|: 

£ = (61,62,63) = (nYi„,{9,^),{rV)Yi^{9,^),{rx V)Yi„,{9,^)) (40) 
For latter reference we mention the following identities: 

V . f = ^(2, -l{l + 1), 0); V X f = i(-e3, 63, -62 - l{l + l)ei) 
r r 

The unperturbed flow is rotation free and the leading contribution to Maxwell's equations [eq. ( p9| 
comes from the perturbed flow 6v = SvnCi + dv±e2, 



(41) 



dr 



dr 
r 



lil + l)B2r = 
l{l + l)E2r^0 



Bi - l{l + l)E3 = 



-D2 7\ = U 

c or 



'-Bs-E, + '-^ = 
c or 



-Ei+lil + l)B3 = ^jl{l + l)f{r, t) 
c 

c ^ dr ^ dr 



r ■ d(rB2) 
-E3 -Bi + = 

c or 



(42) 

(43) 
(44) 

(45) 

(46) 
(47) 

(48) 

(49) 



d6v\ 1 5v\ —6v 



where fir, t) = ^ + 



Notice that B2 , Bi and £'3 are independent of the source term, and are 

(50) 



taken to vanish identically. The other mode is 

E — EiBi + £^262 ; B — B^ez 
Averaging the Poynting vector 



^ c ^ ^ cB* r ^ ^ ^ 

5 = —{E xB*) = —^ -rEi[YVY*) +r'^E2{VY -VY*)!! 



over all directions gives the radial energy flux 



The corresponding spectral intensity is 



cr^l{l + l) 



^2(^)^3(^)1 



(51) 



(52) 



(53) 



We obtain the wave equation for A = Ei{Lu)r by combining eqs.(p^-(p^) 

(V^ + k^)A{uj) ^ ikrjl{l + 1)/(lj, r) 



(54) 
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and in terms of A, the spectral intensity reads 



21(1 + 1) 



rA{uj) 



dirA*{uj)) 



dr 



(55) 



The wave equation is solved through the Green's function method in the region away from the near 
zone: 



A{uj) ^ ikhf-\kr) / [-ikT]l{l + l)f{uj,rl)]ji{krl)r''^drl. 



In the radiation zone, Kiuj) reduces to : 



k{Lo)rK-e'^''{-iy+^k7^l{l + l) I J {uj,rl)ji{krf)rl' dri 



Putting these pieces together, 



Ii{Lo)^\afkH{l + l)\Mk)\^ 



with 



Ai{k)^ J j f{r,t)e-'^']i{krydrdt 
The function f{r,t) can be expressed in terms of the fluctuation functions [eqs.(|l^)]. 



a/3-1 



(56) 
(57) 

(58) 
(59) 

(60) 



Calling X = kr and performing a change of integration variables we obtain radiation emission rate per 
wave-length A: 

(61) 



Pi{\)^^e'a'l{l + l)\Wi{k)\ 
A 



with 



Wi{k) 



poo pi 

/ ji{x) x^dx [^t' + 2T-{l~V)<i>]y"^~^exp[-tQy]dy, 
Jo Jo 



(62) 



where Q = fci?o(^o/-Roc) and p = 0^7^/2 . According to Barber (|^) the ratio aRo/to = cq, the speed of 
sound, and Q = a fci?o(co/c) ~ 10~^(fci?o)- The asymptotic behavior given by eq.(|28|) and the fact that 
V cc ^^1/" suggests the expansion: 



Kr' + 2r - (1 - V)^] = J2 fcn^"^"/" = E ^" 



n=l 



(— 

\kRo 



/3— n/a 



y 



,71 — 0.(3 



(63) 



where the coefficients 6„ are determined by the dynamics of perturbations. Note that the sum does 
not contain the n = term because the leading term of the series [see again eq.(|8| )] for the velocity 
components $, r is 1^"'^""'^. Therefore, 



Wi{k) - UkRoT'"-^ / y"-^exp[-zQy]dy 



n=l 



kR„ 



Jl 



{x) x2+'^-"/-"da;. 



(64) 



The cutoff fei?max in the x-integral was introduced because the shock does not extend beyond i?„ 
the ambient radius of the bubble. For Q << 1 we might transform this expression into 



wm^ikR^r^Y. 



bn 



n=l 



{kRo) 



i/a 



kR„ 
kRa 



Jl 



(x) x^+^-'''^-°'dx + 



kR„ 



' ji{x) x^+'^dx 



(65) 



The detailed form of the spectrum requires a full knowledge of bn, that is to say, dynamics of the fluctu- 
ations must be specified (this can be done analytically by using the method developed in the appendix). 
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Fortunately, the major features of the spectrum can be obtained without delving into the differential 
equations. For instance, in the region where fci?max < 1 we can approximate ji{x) ~ {2xyi\/{2l + 1)1 and 
then 



Wi{k) 



2H\ 



{2l + l)\ 



{kRc 



1 + 3 + 13- n/a 



Rq 



l+3+p-n/a 



il + 3 + P) \ Ro 



Rn 



1+13+3' 



(66) 

In the other end of the spectrum kRo > 1, taking the asymptotic expression ji{x) ~ 1/ x sin{x — Itt /2) 
is justified, either because in the first integral the integration variable x > 1 or because in the second 
integral the measure x'^^^ (with 2 + f3 > 1) ensures that important contributions to the integral comes 
from the large arguments. Thus, 



Wi{k) ^ (kRo) 
where 



R„ 



Rf) 



/3+1 



/(/3; fci?max) + 



Ro 



n/a 



f{(3 - n/a] fciJmax) - /(/3 - n/a; kRo) 

(67) 



fiP;x) 



Im 



-ilTT/2 



E 



{ix) 



m+1 



{m + [3 + 2 — n/a)m\ 



m=0 



(68) 



The dominant power low contribution to Wi{k) in the region kRo > 1 comes from the linear term 
(kRo) because the series f{(3;x) behaves nearly like sin(a::), for x > 1. Taking the following figures 
-Rmax ^ 5/im , the ambient radius of the bubble and Ro ^ 0.15fim, (we shall explain in a moment) and 
defining Aq — 27ri?o, we display our asymptotic expressions in the form 



where 

Ai = l{l+1) 

and 



2H\ 



1)! E 



(2l + iy.^n ] 1 + 3 + P -n/a 
gi(X) = l{l + \) 



Ro 



l+3+P-n/a 



1 



{1 + 3 + (3) \ Ro 



Rn 



Rn 



Ro 



/3+1 



hlik) 



with 



h,ik) = Y.- 

n— 1 



/(/?; kR^^^) + 



Ro 



Rn 



nj a 



f{f3 - n/a; fciimax) - /(/? - n/a; kRo) 



(69) 

1+13+3 

(70)' 
(71) 

(72) 



The apparent divergence of gi{X) at large angular momenta [see eq. (|7lJ)] seems to endanger the 
present results. This worry is removed studying the asymptotic behavior gi{X), bearing in mind that in 
this limit /3 ~ (7 — l)/(7 + 1), [?]. This yields that f;i(A) ^ as / — > 00, regardless of the specific 

form of the dynamical coefficients 6„ may take. 



3 Assessment of the Results 

The present SL mechanism relies on very basic pieces of physics, the existence of corrugation instabilities 
in spherical shocks, whose existence is well known. Maxwell's equations and the inertial polarization 



10 



paradigm. As we had the opportunity to explain, this paradigm stems from very elementary physics and 
it has remained hitherto unnoticed only because huge accelerations are required for sizeable polarizations. 
The detection of shock polarization in non-polar liquids would lend an undisputable status to the inertial 
polarization principle . In the transduction of sound into radiation , the flash of light must be coincident 
with a burst of sound since the emission of radiation is caused by corrugation instabilities . According to 
eg. (|50|) , only one field-mode is related to the sonoluminescent light. This mode has a longitudinal electric 
field component Ei, and some experiment must be devised to detect it .The transversal component 
E2 points into the direction of the vector 



and this (weird) polarization should be observed in sonoluminescent light. 

Physics is seldom controlled by cut-off parameters, and we expect the cut-off parameter Rmax ( the 
bubble's ambient radius) to play a marginal role in delimiting the frequency band where light is emitted. 
The main features of the spectrum should be controlled by the remaining parameters: i?o, the radius of 
the shock- wave when it is first formed and the perturbation amplitude e . Thus, i?o should characterize 
the typical wave-length of the emitted light \ k, \q — 2ttRq. Our asymptotic results [eq.(|69|)] confirms this 
feeling. Numerical and theoretical studies of the dynamics of imploding shocks support the picture that 
the bubble collapses at the speed of sound by the time it passes through its ambient radius as the right 
criterion both for shock formation and the existence of SL (jl^- |l^). According to these investigations, 
at lOOps before the bubble reaches its minimum size, a shock wave of initial radius Rq = 0.15fim develops: 
by this time the interface is imploding with 4 to 5 times the ambient speed of sound. With these figures, 
we predict the emitted light to lie in A « Aq = 900nm spectral region, regardless the kind of gas present 
in the bubble; in SL experiments light is observed in the 200nm < A < 800nm interval. According 
to this result, it is legitimate to infer the spectrum in this wave-length interval through the asymptotic 
formula for A ^ Ao[see eq.(^9|)]. How does the particular kind of gas present in the bubble impact on the 
the emission power? The dependence of the emitted light upon the particular type of gas present in the 
bubble stems from two different factors: 

i.) different values of the adiabatic index 7 leads to a different shock-wave and corrugation instability 
dynamics; ii.) different gases have different dielectric permeability e. 

The dielectric nature of the gas is implemented through the replacement E ^ D in the Poynting 
vector, which corresponds to the replacement of |VF/(fc)|^ by ?(fc) |Ty;(A:)|^ , or gi{X)^ e{k)gi{X) . Different 

adiabatic indexes would cause hi{k) to change because both the spectrum of (3 and the dynamical 
coefficients 6„ , depend upon 7. These two conditions will cause a change on the shape of the function gi (A) 

. Assuming that after taking these corrections into account, the function gi{X) still remains marginally 
dependent upon the wave-length (non power law), the overall change produced by different gases in the 
shape on the logartithmic representation of the spectrum InP ^ — 31nA -I- ln g((A) -t- const for A < Aq 
, is a displacement of the nearly parallel lines of inclination m = —3. This behaviour is changed as 
we approach the A << Ao region because then the dielectric constant being governed by the plasma 
frequency of the gas, causes the function higi{X) to strongly depend upon A . 

Infering the uncorrected spectra for transmission by the surrounding medium observed by Hillcr in SL 
experiments for bubbles trapping pure noble gases bubbles at 0*^(7 ([^) we infered in = —2.7. For pure 
He, m = —2.5. Inspection of the spectra shows the nearly linear dependence for Ar, He, He^, and Ne. 
The agreement is less accurate for Xe and Kr, for reasons which are presently unclear: it might well be 
the that heavier noble gases cannot be handled with the naive classical Inertial Polarization picture, they 
have too much internal structure and must be handled with a full quantum mechanical approach. The 
spectrum for a mixture of 1% of He and N2 closely resemble the behavior of pure He2(||l^). Differences 
might be credited to the superimposition of the Bremsstrahlung spectrum of free electrons of the ionized 




21 + 1 {l-m)\ 
Att (l + m)l 



e™*^ sin (6*) {imPl" {cos{e))'e ^ - P{ 
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N2 gas in the mixture to the original spectrum, or even the effect of the Inertial-Polarization fields upon 
these electrons. 



Regarding now the intensity of the outgoing radiation, it is governed by the product pe^ . A small p 
would require large corrugation instabilities, invalidating the linear regime approximations. Surprisingly, 
p = Ep/{2Ha) ±2 1.47 x 10"'^^ Watt , imploding shocks are fantastic power stations ! Actually, we have 
to worry to have sufficiently small perturbations to fit the experimental data! Typical power emissions 
are of of the order of 10~^^Watt/nm in the Aq region |l^, calling for an amplitude e ^ 10"^^ or 
Sr = sRq ~ 10~'^^m, which being much smaller then the nuclear dimensions can have only a quantum 
mechanical origin. Now, the radius of the shock at the moment it is formed Rq is governed by the radius of 
the bubble wall Kb, by the time it is collapsing at 4-5 times the ambient speed of sound. The dependence 

of the former on the latter is linear. In a semi-classical approach, it is to be expected that the fluctuations 
on the shape of the imploding shock are also governed by bubble wall fluctuations, e = SRo/Rq = 6Rb/Rb- 
The fluctuations of the bubble interface should be of the order of the bubble's Compton wave-length A;, 
and Sr = {Ro/Rb)Xb ~ Ap/A^, where Ap ~ 10~^^m is the Compton wave-length of the proton and A^ 
is the number of gas atoms trapped inside the bubble , A^ ~ 10^. Thus, in this scenario Sr ~ lO^^^m, 
which is close to the amplitude needed to fit the observed intensity of the radiation. 

One of the most intriguing issues in SL is beyond any doubt the noble gas puzzle: only bubbles 
containing noble gas, even at very small concentrations, glow. What can we say in this respect? Does 
our paradigm shed some light in this direction? Here is a clue. As the bubble collapses and the attains 
supersonic regime the adiabatic heating raises the gas temperature to ~ OAeV (||l^) . The gas is further 
heated when it crosses the shock front, the temperature is increased by a factor Af' . This is more than 
enough to bring diatomic gases to their excited states, but not for noble gases. The dipole contribution 
{'^imn \p\ '^im7i) ' f"^ of the cxcitcd states being much larger than the Inertial Polarization Fields will 
wash away information regarding the latter. A full quantum mechanical calculation should resolve this 
issue. 

There are immense challenges ahead. From the theoretical point of view, one needs to calculate the 
detailed spectrum taking full account of the shock dynamics, study the back reaction of the polarized 
fields upon the dynamics, clarify whether the polarization caused by quantum mechanical transitions of 
a diatomic molecule are the culprits for washing out the Inertial Polarization fields, etc. The immediate 
experimental challenge is to detect Shock Polarization in non-polar fluids. If the effect is confirmed in 
non-polar fluids then it will be very hard to defuse the present transduction mechanism. 



Appendix — A semi-analytical solution of the differential equations for the 
perturbed flow 

In order to solve the set of differential equations we split the matrix into its regular and divergent parts 

A 

iV - K) 

with 



M{V) ^ +BiV) (73) 



A = 





1 


01 


02 


03 


04 


\ 


m{Vc) 
















{dP/dV)v^ 




















\ 

















( 0i(v^) h{y) 03(^) 04(^) \ 

Z 2t/-i Z Z 

3 + /3 -Z(/ + l) 

V -K / 



(74) 



where 0„is a short notation for 0a(Fc) and 0a(F) = 0„(T/)/P(T^) - miy^) j (m{y)(dP I dV)v,(y - Vc)). 
Assuming Biy") and |y(y)) regular functions at the critical point permits the expansions B{y) = 
J2nBniV - Vcy';\Y{V)) = EfeYfe(F - K)^ Substitution into the differential equation yields the 
recurrence formulae: 
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^Yo = (75) 
Y„+i = + B„_„Y,„ (76) 

The matrix A possess three distinct null-eigenvectors: 

(-02,</'l,O,O) 

(-03,0,01,0) (77) 

(-04,0,0,01). 

Associated to each one of these eigenvectors we can construct through the recurrence relations 
\Y^^'> {V)) . The solution of the differential equation is the linear combination |y(V")) = X]i=2 4 l^'*^^)) • 
The fulfillment of the boundary requires that 

\XiV,)) = |y(Fi)) =^[c2Yf) +c3Yf +c4Yf ](Fi -y,)^ (78) 

k 

This equation constitutes a set of four equations for the unknown {ci,(3), which can be solved in a 
perturbational approach in powers (Vi — 14), once the state |X(Vi)} is known. The only missing piece of 
information is the set of boundary conditions for the perturbed fields. 



(2) _ 
) ~ 

(3) _ 
) ~ 

(4) _ 



The boundary conditions for the perturbed flow 

Supersonic motion produces a discontinuity in the fluid flow known as a shock wave or simply shock. Let 
us call V2 and p2 the fluid velocity and density, and C2 the speed of sound behind the shock, as measured 
in the laboratory frame(likewise, the subscript 1 refers to the same quantities in the front of the shock). 
The normal to the shock is n and its velocity in the lab frame is v shock- The discontinuities have to fulfill 
the following conditions at the shock surface p 



n X \vi - Vshock\ 



n-[V2~ Vshockl 

n-[vi - Vshock] 



nX \V2~ Vshock 



P2 



7-1 
7+1 



7-1 
7 + 1 



27 . ^2' 
[n ■ Vshock) 



7 + 1 



(79) 
(80) 

(81) 



In the perturbed fiow the shock front is displaced from 'Sq : r — R{t) = to (So + SY^) '■ r — R{t) — 
5r(t, 6, 0) = 0. The corresponding perturbed normal is 5n — —SJSr, the perturbed shock velocity is 5r 
while the location of the shock itself in self-similar coordinate is 1 + (5^, = Sr / R{t). Accordingly, 



Si 

6r 
Sri 
Svs 



ei?(t)f-j r/,„(^,0) 

\R{t)v)Yi„,{e,cj,) 



ea{l + f3)R{t) ( — 



a/3-1 



Y,r 



(82) 
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The first order corrections to the boundary conditions [eqs. ([79|)-(^ ] are: 



S£_ + Svn — Svs I n + Sv± 



dv 



S(, + 5vn — Svs ] fi + 5v± 



+ 6n X n{v — Vg) 
+ 6n ■ n{v — Vs) 



5c| 



5cl + '^Si 



Inserting eqs.(p^-(p^) into these boundary conditions, yields 

l3Vi - V[ 



$1 

Ai 
5Zi 



l-Vi 

= -Vi ; 

Gi 



-n X bvs — Sn X nVs 



= 



7-1 



{5vs + n ■ SnVg) 



7+1 

2Z(1) Vs (Sn ■ vl + fi ■ 6vs) 



or, equivalently 



\XiV,))^e 



-(/3 + 3)yi/(l-T/i) 
V 2(/3+(l/a-yi)/(l-Fi))/7 / 



(83) 



(84) 



(85) 



Numerical Procedure 

Our procedure for resolving the spectrum of /3 consists of first fitting the unperturbed flow (Z{V), dV{V)) 
by a polynomial in V, from which we extract the matrices A and B{V) as power series in V. Then 
through the recurrence formulae (??) we obtain the expansion coefficients Y„(/3) up to a given order 
and insert then into eq. ([78|) , in conjunction with the above boundary condition |X(Vi)) [eq.(|5|) ] . This 
procedure yields a polynomial equation for /? , which is solved numerically. We display the results for 
7 = 7/5,^=1,2,3,4. 
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